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ABSTRACT 


As part of a continuing effort under Contract HAS 7-^58^ an analytical 
study of catalyzed hydrazine decomposition reaction chamhers was performed to 
assess the effects on the steady-state hehavior of the system of nonunifom 
radial injection and of catalyst hed configurations exhihitlng hoth radial 
and axial nonuniformities . Radial variations in mass flow rate or hed packing 
cause radial temperature and concentration gradients which lead to turbulent 
diffusion of heat and mass in the reactor system. A computer program was 
developed to calculate temperature and reactant concentration dlstrlhutlons 
as functions of axial and radial position in typical hydrazine reaction cham- 
bers. The program is based upon a model of the reactor system which includes 
treatment of the turbulent diffusion of heat and mass in the fre'e-gas phase 
along with heat and mass diffusion within the catalyst particles and between 
the particles and the free-gas phase. Both thermal and catalytic decomposi- 
tion of the reactants are considered:. Included in this report are descrip- 
tions of the turbulent diffusion phenomena^ the reactor model which treats 
these phenomena and typical results calculated on the basis of this model. 
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Analytical Study of Catalytic Reactors 


for Hydrazine Decomposition 
Second Annual Progress Report 
April 1^, 1967 - April l 4 , 1968 

Contract Ho. HAS 7-^58 


SUMMARY 


The Research Laboratories of United Aircraft Corporation ijnder Contract 
HAS 7-458 with the Hatlonal Aeronautics and Space Administration are performing 
an analytical study of catalytic reactors for hydrazine decomposition. This 
second annual technical report summarizes work performed under this continuing 
contract from April 15 , I967 to April l 4 , 1968 . Work during this period has 
included the development of a computer program representing a two-dimensional 
steady-state model of a distributed-feed catalyzed hydrazine decomposition 
reaction chamber. This program was developed to assess the effects on the 
steady-state behavior of the system of nonuniform radial injection and of cata- 
lyst bed eonfiguratlons exhibiting both radial and axial nonimlf ormities . The 
program Is based upon a model of the reactor system which includes a treatment 
of the turbulent diffusion of heat and mass in the free-gas phase along with 
heat and mass diffusion within the catalyst particles and between the particles 
and the free-gas phase. Both thermal and catalytic decomposition of the reac- 
tants are considered. Calculations have been made of temperature and reactant 
concent ration distributions as functions of axial and radial position in typi- 
cal hydrazine reaction chambers for a number of injection profiles and catalyst 
bed configurations . 

The integral equation method, developed during the first year of contract 
effort to describe the diffuslonal processes within the pores of the catalyst 
particles, has been modified to include the effects of heat and mass transfer 
from the bulk fluid, through a stagnant film surrounding the catalyst particles, 
and to the outside surface of the particles . An iterative procedure has been 
developed to solve the implicit integral equations describing reactant concen- 
tration and temperature profiles in the porous catalyst particles. A computer 
program representing this iterative procedure has been used as a subroutine in 
the two-dimensional steady-state program. This subroutine has also been used 
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in the one-dimensional program to define more precisely the kinetics of ammonia 
dissociation hy comparing temperature and concentration profiles calculated using 
this program with availahle experimental information. 

Empirical predictions have been developed of one -dimensional steady-state 
temperature and fractional ammonia dissociation profiles in hydrazine reactors 
packed with Shell 405 catalyst particles. The empirical correlations were 
developed on the basis of many runs made with the steady-state computer program 
developed during the first year of effort on the present contract. It was found 
that fractional ammonia dissociation and bulk fluid temperature are easily pre- 
dicted for a broad range of operating conditions for cases in which most of the 
hydrazine decomposition occurs in the first few tenths of an inch of the reactor; 
this rapid hydrazine decomposition rate is associated with reactors packed with 
particles 25 mesh or smaller for approximately 0.2 inch of reactor bed length. 
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INTBDDUCTIOW 


Under Contract NAS 7-^58, the Eesearch Laboratories of United Aircraft 
Corporation are perfomiing analytical studies to characterize the behavior of 
distributed-feed catalytic reactors for hydrazine decomposition. The specific 
objectives of this program are (a) to develop computer programs for predicting 
the temperature and concentration distributions in monopropellant hydrazine 
catalytic reactors in ■which hydrazine can be injected at arbitrary axial and 
radial locations in the reaction chamber and (b) to perform calculations using 
these computer programs to demonstrate the effects of various system parameters 
on the performance of the reactor. 

Progress previously reported in the first annual report (Ref. l) has 
included the development of computer programs •which describe the steady-state 
and transient behavior of a reactor system in ■which complete radial mixing in 
the free-gas (or liquid) phase ■was assumed. These programs had been used to 
calculate temperature and reactant concentration distributions as functions of 
initial bed temperature, feed temperature, chamber pressure, mass flow rate, 
catalyst size distribution, and axial injector locations. 

During the present reporting period attention has been focused on extending 
the steady-state model to include radial as well as axial variations in temper- 
ature and concentrations in order to permit an analysis of various injection 
schemes and catalyst bed configurations which exhibit radial nonuniformities. 

In addition, a refinement has been made of the method used in the one- and two- 
dimensional steady-state models to describe the effects of diffusional processes 
on reaction rates in porous catalyst particles and sample calculations using 
this modified film and pore diffusion model have been made. Effort during the 
second annual reporting period is described in detail in . succeeding sections 
of this report. 
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Discussion 


The analysis of a hydrazine engine reaction system carried out to date 
pertains to a reaction chamber of arbitrary cross section packed with catalyst 
particles into which liquid hydrazine is injected at arbitrarily selected loca- 
tions . Catalyst particles are represented as "equivalent" spheres with a dia- 
meter taken as a function of the particle size and shape. Both thermal and 
catalytic vapor phase decomposition of hydrazine and ammonia are considered in 
developing equations describing the concentration distributions of these reac- 
tants. Diffusion of reactants from the free-gas phase to the outside surface 
of the catalyst pellets is taken into account. Since the catalyst material is 
impregnated on the interior and exterior surfaces of porous particles, the 
diffusion of reactants into the porous structure must also be considered. In 
addition, the conduction of heat within the porous particles must be taken into 
account since the decomposition reactions are accompanied by the evolution or 
absorption of heat. 

Included in succeeding sections of this report are detailed descriptions 
of (a) the development of an integral equation method describing the effects of 
film and pore diffusion of heat and mass on reaction rates in porous catalyst 
particles, (b) the use of a computer program representing this integral equation 
method as a subroutine in the one-dimensional steady-state model of the reactor 
system to define more precisely the kinetics of ammonia dissociation by comparing 
temperature and concentration profiles calculated using this program with avail- 
able experimental information, (c) the use of the one -dimensional steady-state 
program to develop empirical correlations to predict axial temperature and frac- 
tional ammonia dissociation profiles in hydrazine reactors packed with Shell 
405 catalyst particles, (d) the development of computer program representing 
the two-dimensional steady-state model of the reactor system, and (e) the use 
of the two-dimensional program to calculate the effects on steady-state temper- 
ature and reactant concentration distributions of nonuniform radial Injection 
and of catalyst bed configurations exhibiting both radial and axial nonuniformities. 

Integral Analysis of Processes Occurring within Catalyst Particles 


Catalytic reaction of hydrazine on Shell 405 catalyst particles is so fast 
(Ref. l) that, even at low temperatures, the rate of decomposition of hydrazine 
vapor is controlled by the rate of diffusion of hydrazine from the bulk vapor 
through a stagnant gas film surrounding the catalyst particles to the outside 
surface of the particles. In the case of ammonia, however, film diffusion is 
rapid relative to the rate of dissociation of ammonia within the particles. The 
concentration of ammonia at the surface of the catalyst particles, (cp)g®®3j • 
is therefore fairly close to the ammonia concentrations in the bulk vapor phase. 

Cl ^3* “fhe surface concentration can be calculated, along with the concentration 
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profile in the porous particles, at any eixial location in the reaction chamber 
by solving simultaneously the equations representing film and pore diffusion of 
heat and mass. In describing the diffusion of . mass within a porous pellet, it 
is assumed that changes in the mass density of fluid within the particle are 
negligible relative to changes in concentration of the reacting species. In 
addition, pressure changes within the particle resulting from nonequimolar dif- 
fusion are neglected, as is heat transported by pore diffusion of mass. Assuming 
constant diffusion coefficients. Dp, and thermal conductivities, Kp, the equa- 
tions describing heat and mass transfer within a catalyst particle may be written 
as 


NH 

Dp 


3 


2 NH 
C 


NH3 

'"het 


0 


( 1 ) 


Kp Tp 



NH, 

het 


0 


( 2 ) 


The boundary conditions which consider diffiision of heat and mass through a 
film surrounding a spherical particle are 



where (k^ e-_) represents the rate of diffusion of hydrazine to the particle 

surface. It should be noted here that (Cp)g^2ll4 is approximately zero, reflec- 
ting the fact that the catalytic decomposition of hydrazine vapor is quite rapid. 

As in the particle diffusion model developed in Ref. 1 which neglected 
film diffusion (assuming equal surface and interstitial temperatures and ammonia 
concentrations), the temperature and ammonia concentration within a catalyst 
particle are related by (Ref. 2) 
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(Tp)c 


H 


NH, NH, 


Kc 




NH 3 

'S 


NHj 

Cp 


( 5 ) 


The use of "this relationship enables the reaction rate^ ^het^^-’ written 
as a function of concentration alone Instead of concentration and temperature. 
In this case, however, the reaction rate is a function of two parameters, (T„) 


and (Cp)s 3^ which are yet to he determined. Equation (l) can he solved for the 
concentration at any point in the porous particle in terms of the reaction rate^ 
%et™3^ and the interstitial concentration^ C The solution is derived in 
Appendix I as an implicit integral equation given hy 


p^s 


Gp ( X j 


NH, 


NH:, 


I nh, ^ 
akn “D, 


NH, 


-r 


akf 


NH, 




NH, 

Dp 


/' 


NH, 


■2 


(^p) 


DpN«3 




a k 


NH, 


''het fcp) 


(6) 


Dp 


NH, 


dC 


In order to determine the particle ammonia concentration profile directly in 
terms of the interstitial temperature and reactant concentrations it is neces- 
sary to solve Eq.s . ( 3 )^ (?<•), and (6) simultaneously. Numerical methods to accomplish 
this have been developed and programmed for machine computation.* 


In the special case of negligible film resistance to heat and mass trans- 


fer (i.e., (Q^)s = Tj_ and 
species. 


as 




Cp(x) 


■±), Eq. 

1 _ _L 

X a 


(6) can be written, for any reacting 
|hetl£p) 


/ 


/■[ 


■2 Ulfil 


Dc 


(Cp) 


di 


di 


I 

( ~ "a 

It is Eq. (7) which is used to describe the hydrazine concentration profiles 
within the catalyst particles located in the liquid region of the reaction 
chamber. In this liquid region it is assumed that liquid hydrazine wets the 
outside surface of the catalyst particles so that (cp)g E 2 H 4 _ where 


(T) 


Ci ^2^4 is the vapor concentration in equilibrium with liquid hydrazine at 


temperature 


T- . 


*These methods are described in detail in a computer manual describing to 
potential users the operation of the computer programs representing the one- 
and two-dimensional steady-state models of the reactor system. This manual 
is being prepared as part of the third year of effort under the present con- 
tract. 
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Sample calculations of ammonia concentration and temperature profiles within 
a porous catalyst particle have heen made with and without hydrazine in the inter 
stitial fluid. An illustrative example was considered for the case without hydra 
zine for which the parameters are 


T± 


= 2024 deg R 


Mole Fractions 
in 

Bulk Fluid 


im3 

Hs 

N2 


Dp™3 

H®3 


k™3 


A, 


K. 


P 


= 100 psia 

= 0.19 
= 0.51 
= 0.30 

= 0.35 X 10“^ft^/sec 
= 1405 Btu/lh 
= 5*0 ft/sec 

= 0.29 Btu/ft^-sec-deg E 
= 0.40 X 10"^ Btu/ft-sec-deg R 
= 10-3 ft 


The temperature and ammonia concentration profiles within the catalyst particle 
were calculated by simultaneous solution of Eq.s. ( 3 ), ( 4 ), and ( 6 ) using the 
expression for rj^g^ cited in Ref. 1. These profiles are plotted for this case 
in Rig. 1. The flux at the particle surface, D^{dc ^/c/x) ^ or [cj_ - (cp)g] , 
is easily calculated once the concentration profile is known. For the ammonia 
dissociation case discussed above, the mass flux of ammonia into the catalyst 
particle, normalized by dividing by (k^ Cj^)^3, is plotted as a function of 
bulk fluid temperature in Fig. 2. The reactant concentration profile and then 
the mass flux at the particle surface were obtained for temperatures between 
IJOO and 2J00 deg R; all other parameters were fixed at the same values used 
in computing the profiles shown in Fig. 1. The normalizing factor, (kg 
is the ammonia mass flux which would be obtained if the reaction were controlled 
by the film diffusion of heat and mass. For comparison purposes, normalized 
fl-uxes are also plotted in Fig. 2 for the case in which film resistance to heat 
and mass transfer is negligible (i.e., (cp)g = c^^ and (Tp) g= Tj^) and for the 
case where film and pore diffusion are sufficiently rapid so that the system . 
is controlled by the rate of chemical reaction on the catalytic surfaces (i.e., 
Cp(x) = Cj^ and Tp(x) = Tj^). For the case of negligible film resistance to heat 
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and mass transfer, Eq.. (7) was used to determine the ammonia concentration dis- 
tribution within the particle and then the flux -at the particle surface . 
the chemical reaction-controlled system, the flux was computed as (a/3)rhet 
As shown in Fig. 2, these two cases represent low- temperature asymptotes for 
the general case where the effects of film and pore diffusion on catalytic 
reaction rates are considered. At hi^ temperature this flux asymptotically 
approaches the flux which would he obtained if the reaction were controlled by 
the film diffusion of heat and mass. 


Equations (3), (^), and (6) have also been solved simultaneously for the 
temperature and ammonia concentration profiles in a porous catalyst particle 
for a ease in which hydrazine is present in the bulk fluid and the parameters 
are 


T-; 


Mole Fractions 
in 

Bulk Fluid 




W2% 

™3 

N2 


D. 


HHq 
p J 

H™3 

#2% 

HH3 


k, 


k, 


112% 


K. 


= 2000 deg R 
= 100 psia 
= 0.1t3 
= 0.11 
= 0.23 
= 0.23 

= 0.34 X 10"^ft2/sec 

= l404 Btu/lb 
= -1928 Btu/lb 
= 4.4 ft/sec 
= 3*0 ft/sec 

=0.29 Btu/ft^-sec-deg R 
= 0.40 X 10"^ Btu/ft-sec-deg R 
= 10-3 ft 


Both the temperature and the ammonia concentration distributions are plotted 
for this case in Fig. 3* Similar concentration profiles were calculated f or • 
bulk fluid temperatures between 1700 and 2700 deg R; all other parameters were 
fixed at the same values used in computing the profiles shown in Fig. 3* The 
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fl-ux of ammonia at the particle surface -was then calculated, normalized hy 
dividing hy c^)®'3, and plotted as a function of hulk fluid temperature in 
Fig. h. Here again, at high temperature, the flux asymptotically approaches 
the values 'Which "would he obtained if the reaction were controlled hy the film 
diffusion of heat and mass. At low temperattire, in this ease, film resistance 
to heat transfer remains important because of the heat generated at the particle 
surface by the decomposition of hydrazine. 

Use of Integral Analysis in One -Dimensional Steady-State Model 

A computer program representing the integral analysis ■was used as a sub- 
routine in the one -dimensional steady-state program to define more precisely 
the kinetics of ammonia dissociation by comparing temperature and concentration 
profiles calculated using this program with available experimental information. 
The rate of catalytic dissociation of ammonia can be expressed as (Eef . l) 

^-50.000/Tp ib/ft^-sec (8) 

where the concentrations are expressed in Ib/ft^ and Tp is in deg R. A value 
of 0.3 X 10^^ for a® 3 was obtained in Ref. 1 by comparing experimental data 
with calculated steady-state temperature and concentration profiles. However, 
in the model used to calculate these profiles it was assumed that the temper- 
ature at the surface of the catalyst pellet is equal to the interstitial (bulk 
fluid) temperature at the same axial location. Using the more sophisticated 
integral method described above, the same cases were rerun with the one- 
dimensional steady-state program. The use of this method changes the predicted 
steady-state behavior of the reactor most significantly in regions of the 
reactor system where the bulk vapor temperature is low and the heat generated 
at the catalyst particle surfaces by the decomposition of hydrazine res'ults 
in large differences between the bulk vapor tempemture and the particle sur- 
face temperature. High temperatures at the particle surface result in rapid 
ammonia decomposition which leads to low peak temperatures in the bulk vapor 
phase. Apart from lowering the predicted peak temperatures and ammonia con- 
centrations, the use of the modified film and pore diffusion model did not 
result in any appreciable difference between the recalculated temperature and 
mole-fraction profiles and those presented in Ref. 1. The value of a®-3 which^ 
yielded closest agreement with experimental data was 1 x lG^^(lb/ft3)^* (sec)"" . 
The results of calculations using this value are sho'wn in Figs. 5 through l4.* 
For each of these cases the catalyst bed packing was taken to consist of 25-30 
mesh catalyst particles for the first 0.2 in. and l/8 x l/8 in. cylindrical 

*The sensitivity of these results to changes in a®3 is very similar to the 
sensitivity to this rate constant of the results reported in Ref, 1. 
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pellets for the remainder of the hed. Hiis configuration is referred to in 
the figures as the "standard hed configuration."' The reactor operating con- 
ditions of chamber pressure and mass flow rate were varied, as was the feed 
temperature, 1^, although, in each of these cases, was between approx- 
imately 510 and 530 deg R. 

Temperature distributions are plotted in Fig. 5 for a case in which the 
reactor operating conditions were taken as G = 3*12 Ib/ft^-sec and 
p = ^4-79*5 psia. Also shown in Fig. 5 are temperature measurements obtained 
by Eoeket' Research Corporation (Ref. 3) during the cours'e of engine firings 
under the same operating conditions . The calculated mole-fraction profiles 
for this case, together with experimental values of mole fractions, are shown 
in Fig. 6. The results of similar caleulatlons made for G = 2.43 Ib/ft^-sec 
and P = 1042 psia are shown in Figs. 7 and 8; those made for G = 6.29 Ib/ft^-sec 
and P = 9T4 psia are sho'wn in Figs. 9 and 10 j those made for G = 1.52 Ib/ft^-sec 
and P = 217.9 psia are shown in Figs. 11 and 12, and those made for G = I.5I 
Ib/ft^-sec and P = 111.4 psia are shown in Figs. 13 and l4. 

As noted in Ref. 1, the heat transfer coefficients used in this model 
were estimated using correlations developed for nonreacting systems. In regions 
where significant chemical reaction takes place, such as the zones in which 
both hydrazine and ammonia are decomposing, actual heat transfer rates are con- 
siderably higher than those calculated in this model. Precise accounting of 
this increase in rate would lead to a value of a®3 somewhat less than 1 x 10^^, 

The effects of various reactor operating conditions on steady-state tem- 
perature profiles are illustrated more specifically in Figs. 15 through I8. 

A reference case was chosen in which the operating conditions were G = 3*0 
Ib/ft^-sec, P = 100 psia, and Tp = 530 R, and the catalyst bed configuration 
was the "standard bed configuration." Chamber pressure, mass flow rate, 
catalyst bed configuration and axial Injection profile were then varied in 
turn and the resulting temperature distributions were plotted. In Fig. 15 
temperature distributions are plotted for chamber pressures of 100, 500 and 
1000 psia with all other conditions taken as those of the reference case. 
Increasing pressure causes the peak temperature to rise and shift slightly 
toward the inlet of the reactor. This is due to the inhibiting action of 
hydrogen on the rate of ammonia dissociation. Temperature profiles are shown 
in Fig. 16 for mass flow rates of I.5, 3-0 and 6.0 Ib/ft^-sec. Increasing 
flow rate causes the peak temperature to rise and shift slightly away from the 
reactor inlet. The effect of changing the catalyst bed configuration on the 
temperature distribution is shown in Fig, 17* Temperature profiles are plotted 
for beds packed with all 25-30 mesh particles, all 1/8 in. x I/8 in, cylin- 
drical pellets, and the standard nonuniform particle size distribution. It 
is apparent that the larger particles slow down the rates of the catalytic 
decomposition reactions . The effects of distributed injectors on temperature 
profiles are Illustrated in Fig. I8. Temperature distributions are plotted 
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for the reference case, for the case in which 2/3 of the hydrazine is injected 
at the inlet and the remaining 1/3 is injected unifoimly over the first I/2 in. 
of the reactor, and for the case in which 1/3 of the hydrazine is injected at 
the inlet and the remaining 2/3 is injected uniformly over the first I/2 in. 
of the reactor. In these cases the hed was taken to he packed with all 25-30 
mesh particles. 

One -Dimensional Steady-State Program 


A series of runs was made with the one-dimensional steady-state computer 
program in order to develop empirical correlations to predict axial temperature 
and fractional ammonia dissociation profiles in hydrazine reactors packed with 
Shell 405 catalyst particles . Empirical correlations were developed on the 
basis of about 65 nms representing different combinations of mass flow rates, 
pressures and catalyst bed configurations. It was found that fractional ammonia 
dissociation and bulk fluid temperature could be predicted using the equations 

I - Fractional Ammonia Dissociation = <J> (9) 


and 


Tj = 1020 1 $ + [0.075 (P/IOOO)] I + 1535 


( 10 ) 


where 


= (0 66) (G/z)°-^® |[(0.55 0°'’'' - 0.17) (I000/P)°^^] + O.It} (ll) 

and z and a are expressed in ft, G in Ib/ft^-sec, P in psia, and in deg R. 
Here, as in Ref. 1, fractional ammonia dissociation is determined on the basis 
that, with no ammonia dissociation, one mole each of hydrogen and nitrogen are 
formed for every two moles of hydrazine decomposed. Results obtained using these 
equations are illustrated in Figs. 19 and 20, respectively, for cases in which 
most of the hydrazine decomposition occiirs in the first few tenths of an inch 
of the reactorj this rapid hydrazine decomposition rate is associated with 
reactors packed with particles 25 mesh or smaller for approximately 0.2 in. 
of bed length. For these cases the correlations depicted in Figs. 19 and 20 
work well for axial distances greater than one inch and for values of pressure, 
P, between 10 and 1000 psia, mass flow rate, G, between 1.44 and l4.4 Ib/ft^-sec 
(O.Ql and 0.1 Ib/in.^-sec, respectively) and equivalent spherical radius, 
.between 0.001 and 0.01 ft. For a reactor packed with small (< 25 mesh) par- 
ticles for the first few tenths of an inch and larger particles thereafter, 
the particle radius , a , refers to the larger particles . 
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In Figs , 19 20 , locket leseareh e25i)erxmental data (lef . 3) are plotted 

along ■with the empirical predictions and the results of sample cases run using 
the one -dimensional steady-state program. Values of fractional ammonia dis- 
sociation obtained from the steady-state program are plotted for axial locations 
between 1 and 6 inches while values of bulk fluid temperature obtained from the 
program are plotted only for axial locations between 3 and 6 inches. Calcu- 
lated values of temperature for axial locations between 1 and 3 inches scatter 
slightly about the predicted line in Fig. 20 in the same manner as the corres- 
ponding values of fractional ammonia dissociation scatter about the predicted 
line in Fig. 19- 

It should be emphasised that these empirical correlations do not correctly 
predict the behavior of reactors in which hydrazine decomposition is slow, for 
example, reactors which are uniformly packed -with large catalyst particles, such 
as 1/8 in. X 1/8 in. eylinders. Ihe correlations work q.ui'fce well tho-ugh for 
catalyst bed configurations consisting of 25-30 mesh particles for the first 
0.2 in. and 1/8 in. x 1/8 in. cylindrical pellets for the remainder of the bed. 

Two-Dimensional Steady-State Model 

In developing the two-dimensional model of a hydiazlne reactor system the 
temperature and reactant eoncentrations in the bulk fluid phase are permitted 
to vary with radial and axial position in the reaction chamber. In the entrance 
region of the reactor, where the temperature is low enough to permit the exis- 
tence of liquid hydrazine, radial mixing between adjacent layers of liquid is 
neglected. The equations representing the change in liquid enthalpy and tem- 
perature with axial distance at any radial position are the same as those devel- 
oped for the one-dimensional model described in Ref, 1. As in the one -dimensional 
model, catalytic reaction is assumed to be fast enough to keep liquid hydrazine 
from wetting the pores of the particles; the hydrazine concentration at the sur- 
face of the catalyst particles at any location in the entrance region is then 
computed from the vapor pressure of liquid hydi^zine in the interstitial phase 
at the same location. 

In the vapor regions of the reaction chamber, turbulent diffusion of heat 
and mass is considered as a mechanism for radial mixing. Radial heat and mass 
fluxes are computed as functions of temperature and reactant concentration gra- 
dients. Heat is being supplied to the system by homogeneo’us as well as hetero- 
geneous decomposition of hydrazine, and is being removed from the system by the 
catalytic decomposition of ammonia. The change in enthalpy with axial distance 
at any radial location is related to the reactant concentrations in the inter- 
stitial phase and at the surface of the porous catalyst particles by 
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dhj 

di 


|F(hj-hF) + Aphc [Tj-(Tp)s] + 


(I2f 


Ilie changes in reactant ■weight fractions in the interstitial phase ■with axial 
distance at any radial location are related to the reactant concentrations in 
the interstitial phase and at the surface of the porous catalyst particles hy 






I / I NjH 4 ^ Ap , 

^ ■' T \T 'horn ^ ^ 2 


( 15 ) 


+ 


Ap 

2 



Ci ' (Cp)s]) 


NH3 

1^3 




F 



dz 


I 

G 


rj_ 

1 2 M'^2^^4 


(kcCi) 


N 2 H 4 


M ^2^4 


3Ap 

2 


(kc[Ci-(Cp)sj) 


NH3 M 


Hg 


MNH3 


aiMr^^ 

dr 


S 




(16) 


*Equations of this type are presented in some'what different form in Ref. 4. • 
The last term on the right-hand side of the equation reflects the heat 
transferred hy the radial diffusion of mass. 
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where 

Qr - - X OTj /dr ) 


Nr = -€((3ciVdr)^ 


= 0.74 

G 

-0.41 

'Cf g' 




. Pi . 


0.61 G‘ 



-0.e67 

G 

. Pi . 




.ApyU. _ 


-0.41 


The eddy conductivity and diffusivity may he estimated from (Ref. 5] 

OCfG 


and 


€ = 


aG 


58 

The changes in reactant concentrations with axial distance are then 


where 


dCi' 

dz 


Pi 


dWj 

dz 


£l_ ^ 

Pi dz 


£Pi 

dz 


= Pi 


_L 

M dz 


_L ^ 

Tj dz 


I dp 
~P^ 


_L ^ 

M dz 


I 


I dwj'^ 


E ( WjJ/MJ) J ** 
J 
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(19) 

( 20 ) 
( 21 ) 

given hy 

( 22 ) 

(23) 

(24) 
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and the pressure drop may be estimated from the Ergun equation (Ref. i)-) as 

d9 ■ /i-S\ 

dz ^ " V S^j 

The mass flow rate, G, is computed as a function of the rate of feed of liquid 
hydrazine from the distributed injectors into the system. Radial presstire 
gradients, caused by particle-fluid viscous interaction, are neglected. Such 
pressure gradients would lead to recirctilating flow patterns in the reaction 
chamber. This biilk radial flow is neglected in this analysis. It is assumed, 
therefore, that downstream of the injectors the mass flow rate profile remains 
unchanged. 

The equations representing the two-dimensional model have been programmed 
for digital computation . A computer manual containing a detailed description 
of the operating charaeteristies of the program is currently being prepared as 
part of the third year of effort under the present contract. 

Results of Calculations Using the Two-Dimensional Steady-State Program 


A series of calculations ■was made using the two-dimensional steady-state 
computer program in order to examine the effectiveness of the two-dimensional 
model and to evaluate the effects on system performance of nonuniform radial 
injection and of catalyst bed configurations exhibiting both radial and axial 
nonuniformities. The calculated results illustrated in Figs. 21 through 5^ 
refer to a reactor 3 in. in diameter into which liquid hydrazine is injected 
at a temperature of 530 deg R. 

Axial temperature profiles at various radial locations are plotted in 
Fig. 21 for a case in which a radial nonuniformity in mass flow rate, G, is 
represented as a step function (see Fig. 21). In this case the upstream 
chamber pressure ■was taken as 100 psia and the catalyst bed packing ■was taken 
to consist of 25-30 mesh catalyst particles for the first 0.2 in. and l/8 x 
l/8 in. cylindrical pellets for the remainder of the bed. This configuration 
is referred to in the figures as the ’’standard bed configuration". Turbulent 
diffusion of heat, which tends to reduce radial temperature gradients, is more 
prcmoimced in the do^wnstream end of the reactor. Here the ca-fcalyst particle 
size is larger, and both eddy conductivity and eddy diffusivity are directly 
proportional to particle size. The consequences of radial heat conduction 
are complicated somewhat by the simultaneous turbulent diffiision of mass. 
Higher temperatures are associated with more hydrazine decomposition; thus 
high temperature regions may lose heat by mdial conduction, but may gain 
hydrazine from adjoining low temperature regions by radial diffusion of mass. 
Subsequent decomposition of this hydrazine may lead to even higher temper- 


1.75 + 


ISO (I - S ) 
ZaG/fi 


)( 


2aAgc 


( 25 ) 
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atures. In these same regions the ammonia produced hy the decomposition of 
hydrazine may exist at higher concentration than in adjacent low temperature 
regions. In the absence of radial diffusion the ammonia in these high- 
temperature regions would decompose and lower the temperature. With radial 
diffusion, however, the concentration of ammonia available for decomposition 
may be lowered considerably. For the case considered here, these conbined 
effects lead to the temperature distribution shown in Fig. 21*. For comparison 
purposes, the axial temperature profile corresponding to a radially uniform 
mass flow rate of 3-0 Ib/ft^-sec is also plotted in Fig. 21. Biis is the aver- 
age mass -flow rate calculated by averaging the actual mass flow rate profile 
over the cross-sectional area of the reactor. The mole-fraction distributions 
of hydrazine and animonia associated with the temperature distribution shown in 
Fig. 21 are illustrated in Figs. 22 and 23 respectively. 

The results of similar calculations made for two different upstream chamber 
pressures are shown in Figs. 2k- through 29. The calculated temperature, mole 
fraction of hydrazine, and mole fraction of ammonia distributions are plotted 
in Figs. 2h through 26 respectively for a pressure of 200 psia and the same 
mass flow rate profile and bed configuration used in the calculations presented 
in Figs. 21 through 23. These calculations were repeated for a chamber pres- 
sure of 1000 psia j the results are shown in Figs. 27 through 29. 

Retaining the catalyst bed configuration and step function injection pro- 
file discussed above, and taking the upstream chamber pressure as 200 psia, 
similar calculations were made for two other average mass flow rates. The 
calculated temperature, mole fraction of hydrazine, and mole fraction of ammo- 
nia profiles are plotted in Figs . 3Q through 32 respectively for an average 
mass flow rate of 1.0 Ib/ft^-sec and in Figs. 33 through 35 respectively for 
an average mass flow rate of 6.0 Ib/ft^-sec. 

The results of calculations made for two other injection profiles are 
shown in Figs. 36 through 4l. These calculations were made for the "standard 
bed configuration, " an upstream chamber pressure of 100 psia, and an average 
mass flow rate of 3*0 Ib/ft^-sec. The temperature and mole -fraction distri- 
butions associated with the injection profile illustrated in Fig. 36 are 
plotted in Figs. 36 through 38 and those associated with the continuously 
varying injection profile illustrated in Fig. 39 are shown in Figs. 39 through 
hi. 


*The effects of radial diffusion can be more clearly illustrated with cross - 
plots showing radial temperature profiles at various axial locations in the 
reactor. While radial profiles were not plotted for this case, they were 
plotted for another case discussed later in this section. 
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The effects on temperature and reactant concentration distributions of 
two catalyst bed configurations exhibiting both radial and axial nonuniformities 
are illustrated in Figs, k-2 through hj. For both of these configurations the 
mass flow rate was taken as uniform at 3.O Ib/ft^-sec and the upstream chamber 
pressure was taken as 100 psia. The calculated temperature, mole fraction of 
hydrazine, and mole fraction of ammonia distributions corresponding to the bed 
configuration shown in Fig. k-2 are plotted in Figs. k2 through kk respectively. 
Similar calculations corresponding to the bed configuration shown in Fig. ^5 
are plotted in Figs . through . 

The effects of the simultaneous turbulent diffiision of heat and mass on 
temperature and reactant concentration profiles are more clearly indicated for 
a case in which hydrazine injection is uniform across the inlet face of the 
reactor but additional hydrazine is introduced into the reactor through injec- 
tors Imbedded in the catalyst bed. Calculations were made for an injection 
profile of this type, illustrated in Fig. 48, where the chamber pressure was 
taken as 100 psia and the bed configuration as "standard". In this case the 
buried injectors were taken to distribute hydrazine uniformly for 0<r/E<0.7 
over the first 1/2 in. of the reactor, llie calculated temperatures are plotted 
as a function of axial position at various radial positions in Fig. 48. Hydra- 
zine diffusion from low to high temperature regions results in unusually high 
temperatures at the interface between the high and low flow rate regions. This 
results in the formation of a "thermal sheath" which is more clearly illustrated 
in Fig, 49 which is a cross -plot of the results presented in Fig. 48. Here 
temperature is plotted as a function of radial position at various axial loca- 
tions in the reactor. For comparison purposes, the radial temperature profile 
at the exit of a 3 in. bed with a step function (all inlet) injection profile 
is also plotted in Fig. 49. The mole fraction of hydrazine and ammonia' distri- 
butions associated with the temperature profiles shown in Figs. 48 and 49 are 
plotted in Figs. 50 and 51 respectively. 

The effects of various reactor operating conditions on radial temperature 
profiles at the exit of a 3 in. bed are illustrated in Figs. 52 through 54 for 
the standard bed configuration and for hydrazine injection at the reactor inlet 
only. The effect of upstream chamber pressure on exit radial temperature pro- 
file is shown in Fig. 52 for a step function injection profile with an average 
mass flow rate of 3.O Ib/ft^-see. Very little effect is noted over a ten-fold 
pressure range. A marked effect of average mass flow rate on radial temper- 
ature profile is shown in Fig. 53 for the same step function injection profile 
and an upstream chamber pressure of 200 psia. The effect of inlet injection 
profile on radial temperature distribution is illustrated in Fig. 54 for an 
average mass flow rate of 3.O Ib/ft^-sec and an upstream chamber pressure of 
100 psia. 
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It should he noted here that the injection profiles and catalyst bed 
configurations discussed in this section were chosen simply to illustrate the 
two-dimensional effects which can occur in the hydrazine reactor. No attempt 
was made to perfonu an exhavistive study of the effects of radial variations 
in reactor operating and design parameters on the steady-state behavior of the 
reactor system. It is apparent, however, that a one -dimens tonal model of the 
system, based on parameters averaged over the reactor cross-section, is not 
adequate to describe the behavior of a reactor which exhibits significant 
radial variations in injection profile or bed configuration. For these systems, 
it is necessary to use the two-dimensional model to effectively predict, for 
example, the locations of very high temperature regions such as “thermal 
sheaths " . 
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Ap 

Ci 

Cp 

Cp^ 

Cf 

Cp 

Cs 

D-, 

Dp 

F 

9c 

G 

h 

H 

Kp 

M 

M 


Nr 


P 


■ LIST OF SYMBOLS 
Badius of spherical particle, ft 

Total external surface of catalyst particle per unit volume of Led, ft"^ 
Reactant concentration in interstitial fluid, lh/ft3 
■Reactant concentration in gas phase within the porous particle, Ih/ft-^ 
Equals Cp - (Cp)g , lh/ft3 

Specific heat of fluid in the interstitial phase, Btu/lb - deg R 
Average specific heat of fluid in the interstitial phase, Btu/lh - deg R 
Specific heat of catalyst particle, Btu/lh - deg R 

Diffusion coefficient of reactant gas in the interstitial fluid, ft^/sec 

Diffusion coefficient of reactant gas in the porous particle, ft /sec 

Rate of feed of hydrazine from distrlhuted injectors into the system, 
lh/ft3-sec 

Conversion factor, (ihui/lhf ) ft/sec^ 

Mass flow rate, Ih/ft^-sec 
Enthalpy, Btu/lh 

Heat transfer coefficient, Btu/ft^-sec-deg R 

Heat of reaction (negative for exothemic reaction), Btu/lh 

Mass transfer coefficient, ft/sec 

Thermal conductivity of the porous catalyst particle, Btu/ft-sec-deg R 

Bfolecular weight, Ih/lh mole 

Average molecular weight, Ih/lh mole 

Ite,dial mass flux, Ih/ft^-sec 

Chamber pressure, psia 
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Radial heat flux, Btu/ft^-sec 
r Radial distance, ft 

Rate of (heterogeneous ) chemical reaction on the catalyst, surfaces, 

lh/ft3-sec 

r Rate of (homogeneous) chemical reaction in the interstitial phase, 

hom / r, 

lh/ft3-sec 

R Gas constant, equals 10.73 psia - ft3/lh mole - deg R, or, 

Radius of reactor 

t Time, sec 

t* Actual time minus time required, under steady-state conditions, for 

liquid hydrazine to flow from the reactor inlet to the interface between 
the liquid-vapor and vapor regions, sec 

T Temperature, deg R 

Wj Weight fraction of reactant in interstitial phase 

X Radial distance from the center of the spherical catalyst particle, ft 
Z Axial distance, ft 

a Pre exponential factor in rate equation 

S Interparticle void fraction 

€ Eddy dlffuslvity, ft^/sec 

X Eddy conductivity, Btu/ft-sec-deg R 

jJi Viscosity of Interstitial fluid, Ib/ft - sec 

Pj Density of interstitial fluid, lb/ft3 

Pg Bulk density of catalyst particle, Ib/ft^ 

(j) Defined by Eq. ’ (ll) 

Subscripts 

F Refers to feed 

i Refers to interstitial phase 
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P Refers to gas within the porous catalyst particle 

S Refers to surface of catalyst particle 

Superscripts 

J Refers to chemical species 

L Refers to liq,uld at vaporization temperature 
V Refers to vapor at vaporization temperature 
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APPENDIX I 


DERIVA'nON OF HITEGEAL EQUATIONS REPEESENTHIG THE CONCENTRATION 
PEDFILE OF A REACTANT ¥ITHIN A POROUS CATALYST PARTICLE 
SURROUNDED BY A STAGNANT GAS FILM 


In this section equations are developed to describe the steady-state concen- 
tration profile of a reactant within a porous catalyst particle which is surrounded 
by a stagnant gas film. The reactant concentration profile in the porous particle 
at any location can be fomd as a solution to: 


If the catalyst particles are taken to be "equivalent" spheres of radii a , and if 
concentration Cp* is defined such that Cf^ = Cp ~ Cj , Eq.. (l-l) can be written as 




x^ c^x 

( " dx } 


- r. 


het 


- 0 


( 1 - 2 ) 


where X is the radial distance from the center of a sphere. The boundary conditions 
associated with Eq.. (l-2) are 



0 AT X = 0 


d% 


kc ^ 

— Cp AT X = 0 
Dp 


(1-3) 


Eq,uation (l-2) can be rearranged to get 


_d_ 

d% 





<#> ( X , C p* ) 


(1-4) 


The solution to Eq. (l-4) is most easily obtained by converting it into a Fredholm 
integral equation (see Ref. 6 ) of the form 


cj^tx) = 


or / / - 

X U (x) V lx) - u (x) v(x) 


/ 


°G(x.e)-<^lC.Cp*)a'C 


(1-5) 
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•where u ( )^ ) is a solution of 


. 0 

o^x \ / 


(1-6) 


subject to the condition that 


dc 




dQ ^ * 

dx p 


= 0 


J X = 0 


(I-T) 


and v(x) is a solution of 


i («= w) - ° 


(1-8) 


subject to the condition that 


dcp 

''rff 


d\j * 

W ^P 


-J x= a 


(1-9) 


The Green's function, G (x,C) is given hy 


u(0 vlx) FOR 0<C<x 
■' ulx) v(f) FOR x < C ^ a 


(I-IO) 


The function u(x) can he determined hy first integrating Eq.. (l-6) to get 


u = 


T + 


(i-ii) 


Applying Eq.. (I-T) together with the first of boundary conditions (l-3) to 
Eq. (l-ll), it is found that A| = 0 and 


u = 3, 


(1-12) 
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The function v(x) can be determined in a similar manner by first integrating 
Eq.. (l-8) to get 


V = - -f + 02 (1-13) 


and then applying Eq.. (l-9) and. the second of boundary conditions (l-3) to Eq. 
( 1 - 13 ) to get 


V = A; 


0 kc ~ Dp I 

kc X 


(1-14) 


Equations (l-lO), (l-12), and (l-l4) can now be combined to get 

okc“ Dp I J 


G(x,f) = i 


In addition, 


AgB, 


Ag B, 


L a^kc 

0 k^ ~ Dp 
a^kr 


X 

I 

T 


FOR 0 < C < X 
FOR X < C < a 


(1-15) 


x^ [ u (x) v'(x) - u'(x) v(x) ] 


^2® I 


( 1 - 16 ) 


Equations (1-15) and (l-l6) caii now be substituted into Eq. (l-5) to get' 

, X 


Cp*tx) 


+ 


gkc - Dp _ 
a^kc X 

‘Or. _ 

okc ” Dp 


/ 


/ 


or 


C p { X ) - C j 


a2k( 

j_ _ akc-Dp 
X 


o2k, 


/■[ 


ak, 


a^k 


c J 


4> 

r%2 

. J ^ Dp 

Dp 




(I-IT) 


( 1 - 18 ) 


Equation (l-l8) is an implicit integral equation "which can be solved nvmierically 
to determine the concentration at any point in a porous particle in terms of Cj , 
the concentration in the bulk fluid. 
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G9 1046 1-24 fig. 1 

TEMPERATURE AND AMMONIA CONCENTRATION PROFILES WITHIN CATALYST PARTICLE 

NO HYDRAZINE PRESENT IN BULK FLUID 


SEE TEXT FOR VALUES OF PARAMETERS REQUIRED FOR NUMERICAL SOLUTION 
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TEMPERATURE IN INTERSTITIAL PHASE, T - DEG R 
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STEADY-STATE AXIAL TEMPERATURE PROFILE 
USING MODIFIED FILM AND PORE DIFFUSION MODEL 

P = 479.5 P5IA 
G = 3.12 LB/FT2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 



0.05 0.10 0.15 

AXIAL DISTANCE, Z - FT 



MOLE-FRACTIONS OF REACTANTS IN INTERSTITIAL PHASE 
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STEADY-STATE AXIAL PROFILES 
OF MOLE-FRACTIONS OF REACTANTS 
USING MObiFIEb ' FILM AND PORE DIFFUSION MODEL 


P = 479.5 PSIA 
G = 3.12 LB/FT2- SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 
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2 2 


ROCKET RESEARCH DATA 


AXIAL DISTANCE, Z-FT 





TEMPERATURE IN INTERSTITIAL PHASE, T, - DEG R 
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STEADY-STATE AXIAL TEMPERATURE PROFILE 
USING MODIFIED FILM AND PORE DIFFUSION MODEL 

P » 1042 P5IA 
G =2.43 UB/FT^- SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 



AXIAL DISTANCE, Z - FT 



mole-fractions OF REACTANTS IN INTERSTITIAL PHASE 
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FIG. 8 


STEADY-STATE AXIAL PROFILES 
OF MOLE-FRAGTIONS OF REACTANTS 
USING MODIFIED FILM AND PORE DIFFUSION MODEL 

P = 1042 PSIA . 

G = 2.43 LB/FT2- SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




TEMPERATURE IN INTERSTITIAL PHASE, T. - DEG R 
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STEADY-STATE AXIAL TEMPERATURE PROFILE OF I LE 
USING MODIFIED FILM AND PORE DIFFUSION MODEL MODEL 

P =974 PSIA 
G= 6.29 LB/FT SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 



AXIAL DISTANCE, Z - FT 




MOLE-FRACTIONS OF REACTANTS IN INTERSTITIAL PHASE 
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FIG. 10 


STEADY-STATE AXIAL PROFILES 
OF MOLE-FRACTIONS OF REACTANTS 
USING MODIFIED FILM AND PORE DIFFUSION MODEL 


P = 974 PSIA 
G = 6.29 LB/FT2_ sgc 

STANDARD BED CONFIGURATION (SEE TEXT) 




TEMPERATURE IN INTERSTITIAL PHASE, T. - DEG R 
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STEADY-STATE AXIAL TEMPERATURE PROFILE 
USING MODIFIED FILM AND PORE OIFFUSION MODEL 

P = 217.9 PSIA A 
G = 1.52 LB/FT2 - SEC SEC 

STANDARD BED CONFIGURATION (SEE TEXT) TEXT) 



AXIAL DISTANCE, Z - FT 



N INTERSTITIAL PHASE 
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STEADY-STATE AXIAL PROFILES 
OF MOLE-FRACTIONS OF REACTANTS 
USING MODIFIED FILM AND PORE DIFFUSION MODEL 


P = 217.9 PSIA 
G= 1.52 LB/FT2 - SEG 
STANDARD BED CONFIGURATION (SEE TEXT) 






TEMPERATURE IN INTERSTITIAL PHASE, T. - DEG R 
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STEADY-STATE AXIAL TEMPERATURE PROFILE 
USING MODIFIED FILM AND PORE DIFFUSION MODEL 

P= in.4 PSIA 
G = 1.51 LB/FT2 _ SEC 
STANDARD BED CONFIGURATION (SEE TEXT) 



AXIAL DISTANCE, Z - FT 



MOLE-FRACTIONS OF REACTANTS IN INTERSTITIAL PHASE 
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FIG. 14 


STEADY-STATE AXIAL PROFILES 
OF MOLE-FRACTIONS OF REACTANTS 
USING MODIFTEO FILM AND' PORE DIFFUSION MODEL 


P = 111.4 PSIA 
G = 1.51 LB/FT 2- SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 
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ROCKET RESEARCH DATA 
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G 9 1046 1-24 FIG. 

STEADY - STATE AXIAL TEMPERATURE PROFILES 
FOR VARIOUS CHAMBER PRESSURES 
USING MODIFIED FILM AND PORE DIFFUSION MODEL 

G = 3.0 LB/FT 2 - SEC 

STANDARD BED CONFIGURATION (SEE TEXT) 




TEMPERATURE IN INTERSTITIAL PHASE, T. - DEG R 


STEADY - STATE AXIAL TEMPERATURE PROFILES 
FOR VARIOUS MASS FLOW RATES 
USING MODIFIED FILM AND PORE DIFFUSION MODEL 

P = 100 PSIA 

STANDARD BED CONFIGURATION (SEE TEXT) 
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STEADY - STATE AXIAL TEMPERATURE PROFILES 
FOR VARIOUS CATALYTIC BED CONFIGURATIONS 
USING MODIFIED FILM AND PORE DIFFUSION MODEL 

P = 100 PSIA 
G = 3.0 LB/ FT 2 _ SEC 



FIG. 17 





TEMPERATURE IN INTERSTITIAL PHASE, T: - DEG R 


G9 10461 -24 

STEADY - STATE AXIAL TEMPERATURE PROFILES 
FOR VARIOUS HYDRAZINE AXIAL INJECTION PROFILES 
USING MODIFIED FILM AND PORE DIFFUSION MODEL 

P = TOO PSIA 

BED CONFIGURATION: ALL 25-30 MESH GRANULES 







FRACTIONAL AMMONIA DISSCIATION 
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TEMPERATURE IN INTERSTITIAL PHASE, Tj (z) - DEG R 
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MOLE - FRACTION OF AMMONIA IN INTERSTITIAL PHASE 
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FIG. 23 


STEADY - STATE AXIAL PROFILES OF MOLE - FRACTION OF 
AMMONIA AT VARIOUS RADIAL POSITIONS 


P = 100 PSIA 
R =0.125 FT 


AVERAGE G 


0 5 
r/R 


1.0 


STANDARD BED 
CONFIGURATION 
(SEE TEXT) 



AXIAL DISTANCE, z - FT 
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STEADY - STATE AXIAL TEMPERATURE PROFILES 
AT VARIOUS RADIAL POSITIONS 


P = 200 PSIA 
R = 0.125 FT 
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G 3 


AVERAGE G 

_4 I I 1_ 

0,5 

r/R 


jL_i 
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STANDARD BED 
CONFIGURATION 
(SEE TEXT) 



FIG. 24 



MOLE-FRACTION OF HYDRAZINE 
IN INTERSTmAL PHASE 
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FIG. 25 


STEADY -STATE AXIAL PROFILES OF MOLE - FRACTION 
OF HYDRAZINE AT VARIOUS RADIAL POSITIONS 


P » 200 PSIA 
R = 0.125 FT 


- 

STANDARD BED 

AVERAGE G 

CONFIGURATION 

11 1 1 J 1 - J 

. . 1 (SEE TEXT) 
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MOLE - FRACTION OF AMMONIA 
IN INTERSTmAL PHASE 
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FIG. 26 


STEADY -STATE AXIAL PROFILES OF MOLE - FRACTION OF 
AMMONIA AT VARIOUS RADIAL POSITIONS 



AXIAL DISTANCE , i - FT 



TEMPERATURE IN INTERSTITIAL PHASE, 




MOLE - FRACTION OF HYDRAZINE IN INTERSTITIAL PHASE 
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fig. 28 


STEADY -STATE AXIAL PROFILES OF MOLE - FRACTION OF 
HYDRAZINE AT VARIOUS RADIAL POSITIONS 


P = 1000 PSIA 
Rsp 0.125 FT 



STANDARD BED 
CONFIGURATION 
(SEE TEXT) 





MOLE- FRACTION OF AMMONIA 
IN INTERSTITIAL PHASE 
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FIG. 29 


STEADY -STATE AXIAL PROFILES OF MOLE - FRACTION 
OF AMMONIA AT VARIOUS RADIAL 



AXIAL DISTANCE , z - FT 




TEMPERATURE IN INTERSTITIAL PHASE, Tj - DEG R 
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FIG. 30 


STEADY -STATE AXIAL TEMPERATURE PROFILES 
AT VARIOUS RADIAL POSITIONS 

STANDARD BED 
CONFIGURATION 
(SEE TEXT) 

0.5 1.0 

r/R 



AXIAL DISTANCE, z - FT 





MOLE-FRACTION OF HYDRAZINE 
IN INTERSTITIAL PHASE 


G9 10461 -24 


STEADY - STATE AXIAL PROFILES OF MOLE - FRACTION 
OF HYDRAZINE AT VARIOUS RADIAL IPOSITIONS 



AXIAL DISTANCE . z - FT 





TEMPERATURE IN INTERSTITIAL PHASE , T; - DEG R 
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FIG. 33 


STEADY -STATE AXIAL TEMPERATURE PROFILES 
AT VARIOUS RADIAL POSITIONS 


P= 200 PSIA 
R = 0.125 FT 
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MOLE - FRACTION OF AMMONIA IN INTERSTITIAL PHASE 
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FIG. 35 


STEADY - STATE AXIAL PROFILES OF MOLE - FRACTION OF 
AMMONIA AT VARIOUS RADIAL POSITIONS 


P = 200 PSIA 
R = 0.125 FT 
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MOLE - FRACTION OF HYDRAZINE IN INTERSTITIAL PHASE 
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FIG. 37 


STEADY -STATE AXIAL PROFILES OF MOLE - FRACTION 
OF HYDRAZINE AT VARIOUS RADIAL POSITIONS 


P .3 100 PSIA 
R = 0.125 FT 
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MOLE - FRACTION OF AMMONIA IN INTERSTITIAL PHASE 










MOLE - FRACTION OF HYDRAZINE IN INTERSTITIAL PHASE 
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FIG. 40 


STEADY -STATE AXIAL PROFILES OF MOLE - FRACTION 
OF HYDRAZINE AT VARIOUS RADIAL POSITIONS 



AXIAL DISTANCE, z - FT 




MOLE - FRACTION OF AMMONIA IN INTERSTITIAL PHASE 










MOLE - FRACTION OF HYDRAZINE IN INTERSTITIAL PHASE 
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FIG. 43 


STEADY - STATE AXIAL PROFILES OF MOLE - FRACTION 
OF HYDRAZINE AT VARIOUS RADIAL POSITIONS 


P = 100 PSIA R ^ 0.125 FT G = 3.0 LB/FT^ - SEC 


CATALYST BED CONFIGURATION: 
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MOLE-FRACTION OF AMMONIA IN INTERSTITIAL PHASE 
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STEADY - STATE AXIAL PROFILES OF MOLE - FRACTION OF 
AMMONIA AT VARIOUS RADIAL POSITIONS 


R = 0.125 FT G^3.0 LB/FT* -SE( 

CATALYST BED CONFIGURATION: 
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STEADY -STATE AXIAL TEMPERATURE PROFILES AT VARIOUS RADIAL POSITIONS 


P = 100 PSIA e = 3.0 LB/FT ^-SEC 


CATALYST BED CONFIGURATION: 
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MOLE FRACTION OF HYDRAZINE 
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STEADY -STATE AXIAL PROFILES OF MOLE - FRACTION 
OF HYDRAZINE AT VARIOUS RADIAL POSITIONS 


P = TOO PSIA 


R = 0.125 FT G = 3.0 LB/FT - SI 


CATALYST BED CONFIGURATION: 
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MOLE - FRACTION IN INTERSTITIAL PHASE 
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STEADY - STATE AXIAL PROFILES OF MOLE - FRACTION 
OF AMMONIA AT VARIOUS RADIAL POSITIONS 


P= 100 PSIA 


R = 0.125 FT G = 3.0 LB/FT 


CATALYST BED CONFIGURATION: 






TEMPERATURE IN INTERSTITIAL PHASE , T. ^ DEG R 
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STEADY - STATE AXIAL TEMPERATURE PROFILES AT VARIOUS 
RADIAL POSITIONS FOR BURIED INJECTORS 


FIG. 48 















P = 100 PSI 
R = 0.125 F 
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FIG. 53 


EFFECT OF MASS FLOW RATE ON STEADY - STATE RADIAL TEMPERATURE PROFILE 
AT EXIT OF 3 - INCH BED FOR A GIVEN INLET INJECTION PROFILE 

P = 200 PSIA R» 0.125 FT 

STANDARD BED CONFIGURATION (SEE TEXT) 
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FIG. 54 


EFFECT OF INLET INJECTION PROFILE ON STEADY -STATE RADIAL TEMPERATURE 

PROFILE AT EXIT OF 3 - INCH BED 


P = TOO PSIA R = 0.125 FT 

STANDARD BED CONFIGURATION (SEE TEXT) 





